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The interaction of a weakly bound Rydberg electron with 
an electromagnetic half-cycle pulse (HCP) is described with 
the help of a multidimensional semiclassical treatment. This 
approach relates the quantum evolution of the electron to its 
underlying classical dynamics. The method is nonperturba- 
tive and is valid for arbitrary spatial and temporal shapes of 
the applied HCP. On the basis of this approach angle- and 
energy-resolved spectra resulting from the ionization of Ryd- 
berg atoms by HCPs are analyzed. The different types of 
spectra obtainable in the sudden-impact approximation are 
characterized in terms of the appearing semiclassical scatter- 
ing phenomena. Typical modifications of the spectra origi- 
nating from finite pulse effects are discussed. 

PACS numbers: 32.80. Rm, 03.65. Sq 



I. INTRODUCTION 

Recent experimental jj]-|| and theoretical inves- 
tigations have demonstrated that (almost) unipolar, high 
power electromagnetic field pulses are a useful new spec- 
troscopic tool which is particularly well suited for inves- 
tigating the dynamics of weakly bound Rydberg elec- 
trons. In current experiments the duration r of these 
half-cycle pulses (HCPs) extends from the subpicosec- 
ond to the nanosecond regime. Several features distin- 
guish the interaction of Rydberg electrons with this new 
type of electromagnetic radiation from the interaction 
with conventional laser light or microwave fields. First 
of all, and in contrast to the optical case, the HCP in- 
teracts with the Rydberg electron at every point of its 
orbit. Furthermore, due to its unipolar nature the influ- 
ence of the HCP on the electron can often be described in 
relatively simple and intuitive terms, e.g., with the help 
of the so-called sudden-impact approximation where the 
effect of the HCP is modeled as an instantaneous mo- 
mentum change. Finally, the energy transfer between an 
oscillating electromagnetic field and a Rydberg electron 
is always governed by an approximate energy conserva- 
tion whose energy uncertainty is typically small in com- 
parison with the amount of energy transferred. As HCPs 
are almost unipolar such an approximate energy conser- 
vation does not hold and arbitrary amounts of energies 
can be transferred to a charged particle from zero up to 
a maximum amount of the order of Ti/t. 

These characteristic features give rise to interesting 
novel phenomena which have been explored in a number 
of recent investigations. So far these studies have concen- 



trated on total and energy-resolved ionization probabili- 
ties . Theoretical work in this context has focused on 
numerical approaches that often employ the sudden- 
impact approximation and on one-dimensional models 
IIH which try to capture the most significant physical 
effects occuring along the direction of polarization of a 
linearly polarized HCP. 

In this paper a multi-dimensional semiclassical descrip- 
tion of excitation of a weakly bound Rydberg electron by 
a half-cycle pulse is developed. It quantitatively connects 
the quantum evolution of the Rydberg electron to its un- 
derlying classical dynamics. The presented theoretical 
treatment is nonperturbative and it is applicable to ar- 
bitrary spatial and temporal HCP shapes. In particular, 
the approach is not restricted to the sudden-impact ap- 
proximation. Within this theoretical framework measur- 
able transition probabilities are represented as coherent 
sums of probability amplitudes which are associated with 
corresponding classical trajectories of the excited Ryd- 
berg electron. This method is very accurate numerically 
especially for highly excited Rydberg states with large 
principal quantum numbers and for HCPs whose trans- 
ferred energy is small in comparison with typical ground- 
state ionization energies. It is precisely this dynamical 
regime which usually causes severe numerical problems 
in fully quantum mechanical calculations due to the large 
spatial extension of Rydberg states and the presence of 
the Coulomb singularity close to the nucleus. 

On the basis of this theoretical treatment we demon- 
strate that angular distributions of ionization probabili- 
ties contain a wealth of information about the ionization 
dynamics of the Rydberg electron which cannot be ob- 
tained from energy-resolved ionization spectra alone. In 
the sudden-impact approximation, for example, in which 
the pulse duration of the exciting HCP is small in com- 
parison with the classical Kepler period of the Rydberg 
electron energy- and angle-resolved ionization spectra 
are dominated by oscillatory structures and semiclassi- 
cal catastrophes of the rainbow and glory type. These 
structures can be explained in a clear and intuitive way 
in terms of interferences between probability amplitudes 
of a few classical electronic trajectories. The study of 
the influence of finite pulse durations indicates that some 
parts of these oscillatory structures depend strongly on 
the duration of the exciting HCP. A main purpose of our 
subsequent discussion is thus the systematic exploration 
of these latter effects which cannot be described appropri- 
ately by the sudden-impact approximation. First results 
of our studies have also been presented in Ref . . 
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This paper is organized as follows: In Sec. II a gen- 
eral (multidimensional) semiclassical description of exci- 
tation of a weakly bound Rydberg electron by an HCP 
is developed. This general approach is valid for arbitrary 
spatial and temporal dependences of the HCP. Starting 
from this general approach the simplifications arising in 
the sudden-impact approximation are discussed in de- 
tail. In Sec. Ill we first present a characterization of 
the various types of energy- and angle-resolved ioniza- 
tion spectra that can be obtained in the sudden-impact 
approximation. These spectra may be classified accord- 
ing to the appearing semiclassical scattering phenomena. 
Subsequently, it is investigated how a finite duration of 
the exciting HCP influences the behavior of the spectra. 
The comparison with the corresponding results obtained 
in the sudden-impact approximation demonstrate which 
dynamical aspects are particularly sensitive to pulse du- 
ration effects. Finally, in Sec. IV a brief summary and 
conclusions are given. 



II. THEORETICAL DESCRIPTION 

In the first part of this section a nonperturbative semi- 
classical framework is developed that describes the exci- 
tation of a weakly bound Rydberg electron by an HCP 
and that is valid for arbitrary spatial and temporal pulse 
shapes. In the second part this treatment is special- 
ized to the case of short HCPs whose influence on a 
Rydberg electron can be described within the sudden- 
impact approximation. Hartree atomic units are used 
with h = m e = e = 1 (m e and e denote the electronic 
mass and charge). 



A. General multidimensional semiclassical 
description 

Let us consider a typical process in which a weakly 
bound Rydberg electron is excited by a half-cycle pulse. 
The dynamics of the electronic wavefunction \ip) t is 
described by the time-dependent Schrodinger equation 
id\ip)t/dt — H\ip) t with the Hamiltonian 



H = §[-tV, 



-A(x,<)] 2 + y(x). 



(i) 



Outside the core region the effective potential V(x.) ex- 
perienced by the Rydberg electron is of the form V(x) ~ 
— l/|x| whereas inside it is modified due to the presence 
of the other electrons and of the nucleus. The electro- 
magnetic vector potential of the exciting HCP is de- 
noted A(x, i) and fulfills the relations V • A(x, t) = 
and A(x, t — ► — oo) = 0. 

Initially the Rydberg electron is supposed to be pre- 
pared in an energy eigenstate |no^o m o) with principal 
quantum number no > 1 and angular quantum numbers 
lo and mo- Well inside the classically allowed region and 



for low values of the angular momenta, i.e. lo,mo <C uq, 
this state can be approximated semiclassically by || 

{ X \n l a m ) = A ci (x)[e 4S °( x )-"/ 4 + e -«fc(*)4*r/4] (2) 
with the classical eikonal 

Sb(x) - /' ' dr'p(r', eo) - (l a + 1/2)tt + na. (3) 
Jo 

The local radial electronic momentum in the Coulomb 
potential of the ionic core is given by 



p(r',eo) = A/2(eo + l/r') 
with the initial energy 

e = -[2(n -a) 2 }- 1 . 



(4) 



(5) 



While the rapidly oscillating part of the wave function 
(||) is determined by the eikonal So (x) its slowly varying 
amplitude is given by 



A c ;(x) 



i£W(e,$) 



|x|(no-a)3/V2vrp(|x|,eo) 



(6) 



with the spherical harmonic YJ™°(0,<I>). The spherical 
angles of x are denoted 9 and <E>, respectively. Equation 
(||) is valid for electronic distances r — |x| from the nu- 
cleus which are located well outside the core region. All 
electron correlation effects are localized inside this core 
region which has an extension of a few Bohr radii. They 
can be taken into account within the framework of quan- 
tum defect theory p|JlC(]. In the case of an inert ionic 
core which is considered here for the sake of simplicity 
these core effects can be described by a single quantum 
defect a. This quantum defect is approximately energy 
independent sufficiently close to the ionization threshold. 
Generalizations to more complicated situations in which 
electronic core excitations and channel couplings have to 
be taken into account are possible within the framework 
of multi-channel quantum defect theory ]^,[l0| but will 
not be considered in this work. 

In Eq. (||) the initial state \ip)t=o is represented 
by a sum of contributions each of which involves a 
slowly varying amplitude, i.e. A c i(x), and an exponential 
function with a large, imaginary- valued argument, i.e. 
iSo (x) . Within the semiclassical treatment one associates 
with each of these contributions a three-dimensional La- 
grangian manifold in phase space [H 12 1. The special 
form of the initial condition of Eq. (g) implies that this 
Lagrangian manifold has two branches, namely 

L$ = {(x , p ) e R 6 | po - ±V X0 5 (x ) e R 3 }. (7) 

As So(xo) depends only on the radial electronic dis- 
tance r = |x | one finds p = ±|po|x /|x | as long as 
no 3> lo,rriQ. Semiclassically, the solution of the time- 
dependent Schrodinger equation \ip} t is determined by all 
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solutions x±(i;x ,po) and p±(i,x ,po) of the classical 
equations of motion with Hamiltonian 



; [ P -A(x,t)] 2 + y(x) 



(8) 



whose initial conditions are located in this La- 
grangian manifold L^. The probability amplitude 
(x|^) t is determined by all those classical trajectories 
{x^ } (t; x j) , p^, p2 } (t, x^ , p ( j) )} which reach point x 
at time t. 



i 7 Pq_J' 
i.e. Q 



E 



dx A dy A dz 



dxo A dya A dzo 
(i)i 



-1/2N 



a i[Sj(x,i)- w ^/2]^ x ( ) i)|^^ =o 



(9) 



For each contributing trajectory j this probability ampli- 
tude is determined by three characteristic classical prop- 
erties, namely its classical action 

S,(x,t)= [ dt'C(± {j) (t'),^(t'),t'), (10) 
Jo 



the determinant of its Jacobi field fll J 



dx A dy A dz 



dxo A dyo A dzo 



Dct 



d(x,y,z) 



d{xo,yo,z ) 



(11) 



and its Morse index fij. Thereby £(x,x, t) denotes the 
classical Lagrange function associated with the Hamil- 
tonian of Eq. (g). The determinant of the Jacobi field 
characterizes the stability properties of the contributing 
classical trajectory. The Morse index fi is numerically 
equal to the numbers of zeroes of this determinant times 
their multiplicities fl3[| . 

From Eq. (^|) one can evaluate the probability ampli- 
tude of measuring an ionized electron with final momen- 
tum p(f\ With the help of the stationary phase approx- 
imation one obtains the result 



pW)|V)t^oo= lim(27T 



-3/2 



!|e u 'f-p 1J '-x| 



= E P!j cl) e lls ^ pU))+w ^ x o } )-^^/ 2 l (12) 



This ionization amplitude is represented as a sum of all 
trajectories j which start at positions Xq^ with momenta 
Pq and which assume the final momentum p^ long 



after the interaction with the HCP, i.e., at t 
Eq. ( |l2| ) the quantity 



dpi f) A dpi 1 ' A dp. 
dxo A dyo A dzo 



(/) 



-1/2 



oo. In 



(13) 



denotes the classical contribution of trajectory j to the 
angle- and energy-resolved ionization probability. The 
classical action Sj(p) is given by 



Sj (p) = - f dt x« (t) • dp - (t) - x« (t = o) . p (j) (t = o) . 

The phase contribution originating from the initial state 
is determined by 



W i (x) = ±[5 (x)-7T/4] 



(14) 



where one has to choose the plus or minus sign depend- 
ing on whether the initial radial momentum of the Ryd- 
berg electron is positive or negative. The Maslov index 
fij is equal to the number of zeroes of the determinant 

3^ /\d%> Adp \ i — h 

dxo^dyoAdzo j times their multiplicities In 

terms of the ionization amplitude of Eq. ([l2]) the energy- 
and angle-resolved ionization probability is given by 



dP„ 



deU)dQ 



,(/) 



(15) 



with d£l = sin.0 fd<d fdQ f and with 8/ and <&/ denoting 
the spherical angles of the final momentum p^. 

Equations (12) and ( p~5| ) are main results of this sec- 
tion. They yield a complete and numerically accurate 
(compare with Sec. Ill A) semiclassical description of 
the influence of an arbitrary HCP on a weakly bound 
Rydberg electron that is prepared initially in an energy 
eigenstate. These equations are based on Maslov's mul- 
tidimensional generalization of the JWKB method as ap- 
plied to the time-dependent Hamiltonian of Eq. (Q) [jll| . 
It should be mentioned that there exists also a wealth 
of alternative semiclassical methods to approach initial 
value problems. These methods have been pioneered by 
E. Heller [B and W. H. Miller [|l5| and have been gener- 
alized recently in various directions jl6]-[l8| . These types 
of semiclassical approximations are mainly based on the 
dynamics of Gaussian wave packets |l^| and so far they 
have been applied predominantly to problems with ex- 
plicitely time- independent Hamiltonians. 

If the initial state of the Rydberg electron is not an 
energy eigenstate of the form of Eq. (||) almost all of the 
considerations of this section still apply. It is only the 
classical Lagrangian manifold associated with the initial 
state that has to be modified appropriately. In the case 
of a spatially localized electronic Rydberg wave packet, 
for example, which can be described by a quantum state 
of the form 



(x|V) t=0 = S(x)e 



iR(x) 



(16) 



with a large eikonal i?(x) and a slowly varying amplitude 
B (x) the associated Lagrangian manifold is given by 



L = {(x ,Po) G R 6 |Po = V xo i?(x ) e R 3 }. 



(17) 



Contrary to the initial state considered in Eq. (g) this 
Lagrangian manifold has only one branch so that the 
structure of the interferences appearing in Eq. ( |l2| ) is 
expected to be changed significantly. 
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Finally, we would like to mention a useful scaling prop- 
erty of the classical Hamiltonian (||). If spatial varia- 
tions of the pulses are negligible the classical trajectories 
(x(i),p(i)) and (x(t),p(i)) induced by two HCPs A(t) 
and A(i) — jA^t) are related by the scaling relation 



x(t) = 7" 2 x(7 3 *), 
p(t) = 7P(7 3 *)- 



(18) 



These scaling properties are exploited in Sec. Ill A where 
they simplify the classification of the types of spectra 
obtainable in the sudden-ionization approximation. 



B. Semiclassical treatment in the sudden-impact 
approximation 

If the pulse duration t of the exciting HCP is small 
in comparison with the classical Kepler period T c ; = 
27r(no — a) 3 of the initially prepared Rydberg electron, 
the determination of \tp)t can be simplified considerably. 
Neglecting for the sake of clarity any spatial dependence 
the vector potential of the applied HCP can be approxi- 
mated by 



A(t) = -9(t)Ap 
with the characteristic momentum 

Ap = - / cKE(f). 



(19) 



(20) 



Thereby, E(t) = -dA/dt denotes the HCP field strength. 
Thus the solution \ip) t of the time-dependent Schrodinger 
equation with Hamiltonian (jj]) is approximately given by 



-iH A t e iApxl 



(t>0) 



(21) 



with the atomic Hamiltonian Ha = l/2(— iV x ) 2 + V(x). 
Within this sudden-impact approximation the Rydberg 
electron experiences a sudden change of its momentum 
by the amount Ap and subsequently evolves under the 
influence of the atomic Hamiltonian Ha without being 
further affected by the HCP. In this case ionization prob- 
abilities can also be evaluated quantum mechanically by 
means of a partial wave analysis. 

In the case of a hydrogen atom, for example, which is 
prepared initially in the state \no, Iq = mo = 0) the angle- 
and energy-resolved ionization probability is given by the 
partial wave expansion 

HP 1 °° 

-TT^k = r- Ye^PdcosGf) X 



IP ' I 2=0 

drS Wo=0 (r)jK|Ap|r)i^e (/) ;r)- 



(22) 



Thereby Pi and ji are the Legendre polynomial and the 
spherical Bessel function of order I, respectively. The 



regular, energy normalized radial Coulomb wave func- 
tion is denoted Fi(e^;r) and ai is the Coulomb phase 
shift. The radial wave function of the initially prepared 
Rydberg state \uq,Iq — m o = 0) is denoted 5 no ,; =o- In 
Eq. (^2|) it has been assumed that the HCP is linearly 
polarized and that the spherical angles (0/,$/) of the 
final momentum p'-^ are measured with respect to the 
direction of polarization. 

Alternatively, \ip)t can also be determined semiclassi- 
cally with the approach developed in Sec. II A. For this 
purpose one starts from the electronic state which is mod- 
ified by the sudden momentum change due to the HCP. 
Immediately after the application of the HCP this state 
is given by 



) f=+0 = e* A P- A 



lV>>t=- 



(23) 



If the Rydberg electron is prepared initially in the energy 
eigenstate \ip)t=-o = |"o'o m o) of Eq. (|2|), the classical 
Lagrangian manifold which is associated with this state 
is given by 

Lq = {(x ,p ) G R 6 |p = ±V X0 5 (x ) + Ap G R 3 }. 

(24) 

In analogy to Eq. (g) the probability amplitude (x|-0) t is 
determined by all solutions of the classical equations of 
motion with Hamiltonian 



(25) 



whose initial conditions are located in this Lagrangian 
manifold. As this Hamiltonian is explicitly time indepen- 
dent the energy of these trajectories is conserved during 
their time evolution. This implies that outside the core 
region where V(x) w — l/|x| the initial positions Xo on 
this Lagrangian manifold which yield ionizing trajectories 
with final energy e^) are determined by energy conser- 
vation, i.e. 



i[±V Xo So(x ) + Ap] 2 -^=6^. 

^ Xo 



(26) 



Inserting Eq. (|^) and assuming a linearly polarized HCP 
with Ap = Ape z , Ap > 0, yields 



J_ _ (eW - eo - Ap 2 /2) s 
|x | 2Ap 2 cos 2 6 



(27) 



with 0o denoting the polar angle of the initial position 
xo. Equation ( j27| ) determines the initial positions Xo 
of the classical trajectories which ionize with asymptotic 
energy e^. According to Eq. (E7I) there is one value of 
ro = |xo| for each final energy e™ and initial angle 8o- 
At the critical final energy e~|] t = eo + Ap 2 /2 the local 
radial momentum of Eq. (^J) vanishes and the primitive 
semiclassical approximation for the initial state |noZo TO o) 
is no longer applicable. In this special case a proper semi- 
classical description of the initial state can be obtained 
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with the help of uniform semiclassical methods [M • This 
peculiar property suggests that one can at least distin- 
guish two different energy regimes in the excitation of a 
Rydberg electron by a HCP, namely a region of low final 
energies for which < e£} it and a high energy regime in 
which e^) > e^] t . For a fixed value of the final ionization 
energy e^' the initial positions of the classical trajecto- 
ries which ionize into a particular final angle / have to 
be determined by the solutions of the classical equations 
of motion with Hamiltonian (|2^). In the special case of 
a hydrogen atom these relations can be evaluated ana- 
lytically [fH). A comprehensive discussion of energy- and 
angle-resolved ionization probabilities which is based on 
the semiclassical results of this section is presented in 
Sec. III. 



III. ENERGY- AND ANGLE-RESOLVED 
IONIZATION SPECTRA 

In this section energy- and angle-resolved ionization 
probabilities are discussed on the basis of the semiclassi- 
cal approach developed in Sec. II. In the first part charac- 
teristic predictions of the sudden-impact approximation 
are discussed in the semiclassical limit. The second part 
discusses modifications of these ionization spectra that 
arise from a finite duration of the exciting HCP. 



A. Sudden-impact approximation 

This part gives a general overview of the character- 
istic properties of angle- and energy-resolved ionization 
spectra in the sudden-impact approximation for the case 
lo, tuq <C fiQ. As it will be shown below these spectra can 
be described very accurately with the help of semiclassi- 
cal methods. Thus the characteristic qualitative aspects 
of their behavior may be understood by studying the un- 
derlying classical dynamics. 

The classical trajectories pertaining to the parameter 
sets (eo,Ap, e^) and (7 2 eo, jAp, j 2 e^), 7 > 0, are re- 
lated by the scaling transformation of Eq. ([l8|). There- 
fore it is sufficient to examine the classical dynamics as 
a function of the parameters Ap 2 /2|eo| and e^/|eo|, for 
example. The result of a corresponding study is summa- 
rized in Fig. nl In this diagram various regions in the 

ah. 



(Ap 2 /2|eo| 



|eo|)-plane are identified that give rise to 



different types of behavior of the ionization spectra. 

Region (I) pertains to parameters with < < 
e£} it = Ap 2 /2 + e . According to Eq. (27), in this case 



classical trajectories start with negative (positive) radial 
momenta — (+)V Xo S , o(x ) from the positive (negative) z- 
half plane. In order to illustrate the typical behavior of 
the resulting spectra we consider as an example the pa- 
rameter values Ap 2 /2|e | = 6.25, e (/) /ko| = 1-84. Figure 
||(a) depicts the initial positions of the classical trajecto- 
ries in the (x, z)-plane and their respective final angles 



for the initial principal quantum number no — 50 (which 
implies = 10 meV). To show the distribution of the 
scattering angles clearly, in Fig. ||(b) the "deflection func- 
tion" is displayed, i.e., the dependence of the final angle 
Qf on the polar angle Qi of the corresponding initial po- 
sition. From Figs. y(a) and (b) it is apparent that three 
different initial positions pertain to each final emission 
angle 0/. The relative significance of these trajectory 
classes for the ionization signal can be inferred from Fig. 
||(c) where their classical weights P 2 ; = |P c ;| 2 /|j4 c /(x )| 2 
(i.e., the Jacobi determinants) are depicted as a func- 
tion of the final angle f . For small and very large final 
angles the contributions of all three classes are relevant 
whereas in the intermediate region the weight of class 
(II) becomes negligibly small. This decrease is due to 
the fact that class (II) includes trajectories with initial 
positions Xo in the vicinity of x = and the behav- 
ior of such trajectories depends very sensitively on the 
initial conditions. As shown in Fig. ||(d) the existence 
of different contributing trajectories with approximately 
equal weight leads to pronounced interference effects in 
the ionization spectrum. In particular, for small angles 
one observes rapid oscillations which gradually diminish 
as 0/ is increased. This diminution is directly related 
to the peculiar 6 /-dependence of the classical weight of 
trajectory class (II) that was described above. 

An important feature of Figs. |^(a) and (b) is the ex- 
istence of trajectories with Qf = 0°(180°) whose initial 
positions do not lie on the z axis. Consequently, due to 
the rotational symmetry of the system around this axis, 
there are infinitely many classical trajectories that con- 
tribute to the ionization amplitude at these final angles so 
that the classical probabilities |P C /| 2 diverge there. This 
phenomenon is known as semiclassical glory scattering 
p0| . For a proper semiclassical description the primitive 
ionization amplitude of Eq. (|l2|) has to be regularized 
with the help of uniformization methods |po| . Thereby, 
one replaces the contributions of the trajectory classes 
a and b at both sides of the divergence by the uniform 
expression 



e i[(/.+A)/a-/W2-*/4] J n{fa ~ /6) {( go + g b ) x 
M(U - fb)/2] + i(g a - g b )Ji[(f a - fb)/2]} (28) 



with fj = Sj + Wj(x^) and gj = P^ 1 ' (j =a,b). The 
trajectories are labeled in such a way that f a > ft so 
that jla = jib + 1. The Bessel functions of order zero 
and one are denoted Jo(x) and Ji(x). The contribution 
of the third trajectory can be described by the primitive 
semiclassical expression given in Eq. (|l2|). In Fig. ||(d) 
the quantum mechanical ionization spectrum (full curve) 
is compared to the uniform (dashed) and the primitive 
(dotted) semiclassical approximation. All spectra in this 
section are calculated for an initial hydrogenic energy 
eigenstate with Iq = 0. The uniform curve describes the 
quantum mechanical result very well over the whole an- 
gular range whereas the primitive approximation diverges 



(ci) 
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in small regions around 0/ = 0° and 180°. The char- 
acteristic qualitative features of the example described 
above, i.e., the existence of three contributing trajecto- 
ries for all angles 0° < 0/ < 180°, their relative weights 
and the appearance of forward and backward glory effects 
are characteristic for all spectra represented by region (I) 
of Fig. 0} 

The second region of particular importance in Fig. |l| 
is denoted (II). This region is bounded from below and 
above by the curves 

E 2 + 2(1 - 3A)E + SA^/EA - 3A 2 - 2A + 1 = (29) 

and 

E 2 + 2(1 - 3A)E - 8AVEA - 3A 2 - 2A + 1 = (30) 

with E = e^/M and A = Ap 2 /2|e | §§. For A > 1, 
relations (|2^) and ( |30| ) can be approximated as E = 
A — 1/2 and E = 9A — 3/2, respectively. The spectra in 
region (II) are characterized by the appearance of a semi- 
classical rainbow scattering phenomenon and a backward 
glory whereas the forward glory present in the spectra 
of region (I) has disappeared. Again this behavior is 
most clearly illustrated by studying a particular exam- 
ple. To this end, the parameter values Ap 2 /2\e \ = 6.25, 
e (/)/|e | = 7.35 ar e used (with no = 50 in Fig. || which 
implies e^' = 40 meV). In contrast to the previous case 
(and, more generally, as soon as e^' > eo + Ap 2 /2) tra- 
jectories originating from the positive (negative) z-half 
plane now start with positive (negative) radial momen- 
tum. Figure |^(a) depicts the location of initial positions 
in the (x, z)-plane and their corresponding final angles 
whereas Fig. ||(b) shows again the deflection function. 
Apparently, for final angles larger than the "rainbow an- 
gle" r the contributions of three different trajectory 
classes have to be taken into account. At r , however, 
two of these classes coalesce which gives rise to the semi- 
classical rainbow scattering phenomenon |^o| . Again the 
classical weights of the corresponding classical trajecto- 
ries diverge at this angle [cf. Fig. ||(c)].To determine the 
transition amplitude the primitive approximation has to 
be replaced by a uniformized expression. Using the nota- 
tion of expression ( p8|) and defining £ = [3(/i — /2)/4] 2 / 3 
this uniform transition amplitude takes the form 

^[(/ 1+ / 2 )/2-AW2-V4] {(5l +g2) x 

?/*m(-z) + i( 9l - ff2 )r 1/4 Ai'(-e)}. (31) 

In ( |3l| ) Ai(£) and Ai'(£) denote the Airy function and its 
derivative. Quantum mechanically, the rainbow scatter- 
ing phenomenon does not cease abruptly at the rainbow 
angle but extends smoothly into the classically forbidden 
region. This may be taken into account semiclassically 
by extrapolating the quantities fi, gi, and £ of expression 
(|l|) to values 6/ < 6 r . 

Furthermore, Figs. |](a)-(c) show that the forward 
glory has disappeared but that the backward glory is still 



present. The classical weight associated with the differ- 
ent types of trajectories can be inferred from Fig. ||(c). 
Figure |](d) shows that the most prominent difference be- 
tween the spectra of cases (I) and (II) is that the latter 
start for small angles with a smooth non-oscillating part. 
Oscillations set in only around r . 

If, for a fixed value of Ap 2 /2|eo|, the final energy 
e /l e o| is increased the rainbow angle O r grows towards 
the value of 180° and the interference effects in the spec- 
trum are more and more diminished. If the line given 
by Eq. (^) is crossed, rainbow and backward glory scat- 
tering vanish completely. In this region (IV) the spec- 
trum is determined by the contributions of a single tra- 
jectory class and interferences are no longer visible. How- 
ever, Fig. [I] also indicates that in this dynamical region 
the angle-integrated ionization probability dPi on /de^ is 
very small so that this type of behavior is only of little 
practical relevance, in general. 

Regions (I) and (II) are separated by a small area (III) 
in which the classical dynamics is more complicated. In 
this transition region besides the backward glory both 
forward glory and rainbow scattering characteristic of re- 
gions (I) and (II) are present. In addition, however, a sec- 
ond classical rainbow angle appears at a very small value 
of O/. A proper semiclassical description of this behav- 
ior would require uniformization based on a higher-order 
catastrophe p0[ . 

Figures |^(a) and |3|(a) indicate that there is a strong 
correlation between the final direction into which the 
electron is emitted and its position before application 
of the pulse. This correlation might be useful for prob- 
ing the spatial distribution of localized wave packets. In 
order to assess the perspectives of such an approach we 
plot in Fig. |] several ionization spectra obtained in the 
sudden-impact approximation. They pertain to a ra- 
dial Rydberg wave packet detected at different stages of 
its time evolution. For the calculation, the wave packet 
was assumed to be generated by a weak and short laser 
pulse as described in Ref. [^3| and the electronic angu- 
lar momentum was taken as Iq — to show the effects 
in question most clearly. The spectral envelope £(e) of 
the laser pulse was chosen as £ (e) oc exp[— (e — e) 2 (5t) 2 ] 
with the mean excited energy e — —2 x 10~ 4 a.u. corre- 
sponding to a mean quantum number n = 50 and an 
energy spread determined by St — 0.035 T c \. Figure 
0(a) shows the spectra (on a non-logarithmic scale) for 
gU) = 40 meV and various delay times between the ex- 
citing laser pulse and the HCP with Ap — 0.05 a.u. All 
delay times are chosen such that the wave packet has 
not yet reached its outer turning point and is still mov- 
ing outward. The inset depicts some corresponding radial 
distributions IV-'WI 2 f° r the wave packet. The connection 
between these distributions and the form of the spectra 
can be made with the help of the right part of Fig. ||(a) 
where z > 0. As long as I^WI 2 is peaked at a distance r c 
less than ca. 4000 a.u. the spectrum attains its maximum 
at an angle 9/ > whose value agrees with what can be 



G 



inferred directly from r c and Fig. ||(a). For larger values 
of r c the maximum is reached at 0/ = 0° but decreases 
as the wave packet moves beyond 4000 a.u. 

For the case of an inward moving wave packet which 
is examined in Fig. || the spectra are essentially peaked 
around the rainbow angle O r independent of the delay 
time. This observation agrees again with the classical 
picture obtained from Fig. || as the associated weight of 
the electron trajectories diverges at this angle. In or- 
der to achieve a connection between the maxima in the 
spectrum and the electronic probability distribution sim- 
ilar to Fig. ||(a) one should investigate spectra for final 
energies e^) < e£} u [cf. Fig. @(a)]. In summary, Figs. 
H demonstrate that the spatial distribution of the wave 
packet is indeed reflected in the ionization spectrum. In 
order to obtain a complete picture of the distribution 
one will have to analyze a set of spectra with different 
parameters. A particular advantage of the use of HCPs, 
compared, e.g., to optical methods, is that they can in- 
teract with the wave packet at every point of its orbit. 



B. Finite-duration HCPs 

In this subsection it is discussed to which extent the 
results of Sec. Ill A can be applied to the study of angle- 
and energy-resolved ionization spectra from HCPs of fi- 
nite duration. It is natural to examine this problem 
within a semiclassical approach as the theoretical frame- 
work discussed in Sec. II can be applied equally well to 
both instantaneous and time-dependent pulses and its 
quantitative accuracy has been demonstrated in Sec. Ill 
A. Thus the main task consists in studying the classi- 
cal electronic dynamics. In a fully quantum-mechanical 
treatment the partial wave expansion of Sec. II B would 
no longer be applicable. One would have to turn to a 
numerical integration of the time-dependent Schrodinger 
equation which is significantly more complicated due to 
the Coulomb singularity at the nucleus and due to the 
large spatial extension of highly excited Rydberg states. 

In the present context our main question is to which 
extent the essential results of the sudden-impact ap- 
proximation can be recovered in experiments with finite- 
duration HCPs. It is therefore reasonable to concentrate 
the study on short pulses with durations r up to the or- 
der of, say, one tenth of the Rydberg electron's orbit time 
T c i. Typically this corresponds to pulse lengths around 1 
ps which can easily be produced in the laboratory. Before 
discussing in detail an example that illustrates character- 
istic effects originating from finite-duration HCPs let us 
briefly summarize the main results of our investigations: 

(i) The behavior of trajectories starting close to the nu- 
cleus is modified significantly even by very short pulses 
whereas the dynamics of orbits originating at a large dis- 
tance from the core is changed only very gradually. This 
is of course due to the different initial velocities of the 
trajectories: far way from the nucleus the electron moves 



slowly and the pulse can still be considered as almost in- 
stantaneous. From Figs. || and || it follows that the part 
of the spectrum with small final angles 8/, where the 
ionization probability is concentrated and where forward 
glory and rainbow scattering can be observed, is related 
to trajectories with large initial distances from the nu- 
cleus. Thus these essential regions of the spectrum can 
be expected to be relatively insensitive to effects of finite 
pulse durations. On the other hand, the behavior of the 
spectrum at large final angles is determined by trajecto- 
ries starting close to the core. In this region of the spec- 
trum, where the ionization probability is relatively small, 
pulse duration effects are most prominent. In particular, 
in our examples we even find backward glory scattering 
to vanish. 

(ii) The deviations of the spectrum from the sudden- 
impact approximation are mainly determined by the 
pulse duration whereas details of the pulse shape seem 
to be of less importance. We compared the classical dy- 
namics and spectra resulting from two types of HCPs 
the shapes of which were almost rectangular and a half- 
period sine wave, respectively. These choices should rep- 
resent typical variations of realistic HCP shapes. The 
differences in the resulting spectra for these two pulse 
forms (calculated for identical pulse durations and inte- 
grated field strengths) were small in comparison to the 
deviations from the sudden-impact approximation (see 
Figs. |(c) and |). 

(iii) In the sudden-impact approximation there are 
classical trajectories for all values of final energies < 

< oo. With increasing pulse duration, however, we 
find that the range of accessible final energies becomes 
more and more reduced in accordance with the energy- 
time uncertainty principle mentioned in the introduction. 

In order to discuss effects of finite pulse duration in 
more detail we consider as an example the ionization 
spectrum for an HCP with an almost rectangular pulse 
shape given by E(t) = e z i?oexp[— (t/to) 8 ]. The parame- 
ters Eq and to are determined so that the integrated field 
strength and the pulse duration fulfill Ap 2 /2|eo| =6.25 
and t/T c i = 0.05, respectively. Again, for the calcula- 
tion an initial hydrogenic energy eigenstate with Iq = 
and with principal quantum number no — 50 was chosen. 
Figure |5| shows an analysis of the corresponding classical 
dynamics and of the resulting ionization spectrum for a 
final energy of e^'/|e | = 7.35 (compare to Fig. ||). Fig- 
ure ||(a) details the changes in the classical dynamics. 
Three major modifications can be observed: 

(i) The final angles 0/ of trajectories with initially 
outgoing radial momentum do not cover the full range 
0° < 0/ < 90° as in the sudden-impact approximation. 
Instead, they are restricted to values < 0/ < O m with 
6 m w 23.7° in the present case. Here, this does not lead 
to major changes in the spectrum. This is because in 
the angular region where the deviations of this trajectory 
class from the sudden-impact case become significant the 
dominant contributions already originate from trajectory 
classes (I) and (II) [see Figs. ||(b) an d ||(b)]. However, 
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due to the absence of trajectory class (III) for O f > m 
the small short-period modulations of the spectrum of 
Fig. ||(c), which are visible for Qf between 25° and 55°, 
are no longer present in Fig. ||(c). 

(ii) A new class of trajectories appears which has no 
counterpart in the sudden-impact approximation. These 
orbits have initially incoming radial momenta and orig- 
inate from the half-plane z > 0. Their final angles vary 
between 0° and 27.5°. The maximum final angle is not 
attained at [x = 0, z = 0) but at a finite distance from 
the nucleus so that there is a rainbow effect connected 
with this class of trajectories [cf. Fig. ||(b)]. However, 
due to the low classical weight of these trajectories they 
show up in the spectrum most prominently in the form 
of small modulations in the vicinity of 8/ = 0° . 

(iii) The initial positions of trajectories which start 
with incoming radial momenta from the half-plane z < 
are changed significantly. In addition, the final angle 
0/ = 180° is reached only for initial positions with z = 0, 
so that backward glory scattering has disappeared. The 
changes in the ionization spectrum at large angles which 
are due to these effects represent the most significant in- 
fluence of finite pulse durations. 

In Fig. ||(c) the resulting ionization spectrum (full 
curve) is compared to a spectrum which has been 
calculated for an HCP of identical integrated field 
strength and pulse duration but with shape E(t) = 
e z E' sin(27rf/r'), t < r'/2 (dashed curve). The spectrum 
for the sudden-impact approximation is also shown (dot- 
ted curve). As mentioned above, for these pulse durations 
the influence of the HCP shape is only secondary even 
for large final angles. However, as the pulse length is in- 
creased to values larger than ca. 0.1 T or f, these influences 
eventually become more significant. Figure || shows the 
results of an investigation similar to Fig. ||(c) for a final 
energy of e^/kol = 3.68 (e^ = 20 meV) which corre- 
sponds to region (I) of the sudden-impact approximation 
(cf. Fig. |l|). As in the previous case, it is found that 
the main features of the spectrum of the sudden-impact 
approximation are left almost unchanged. The classical 
dynamics change similarly to the previous example. In 
particular, the new trajectory class still exists, but in the 
rapidly oscillating part of the spectrum its contribution 
is almost undiscernible. 



IV. SUMMARY AND CONCLUSIONS 

A general semiclassical treatment of excitation of 
weakly bound Rydberg electrons by intense HCPs has 
been presented. Thus a quantitative connection between 
observable transition probabilities and the underlying 
classical dynamics of the excited Rydberg electron has 
been established. This approach is numerically accurate 
for high principal quantum numbers of the Rydberg elec- 
tron and in cases in which the energy transfer from the 
exciting HCP to the Rydberg electron is small in compar- 



ison with typical ionization energies of low-lying states. 
The underlying classical dynamics of the Rydberg elec- 
tron yields a clear and natural explanation for the os- 
cillatory structures which govern the resulting energy- 
and angle-resolved ionization spectra. These structures 
arise from interference between probability amplitudes 
associated with different ionizing trajectories. Thereby 
two types of semiclassical catastrophes appear, namely 
glory scattering for electron emission in the forward and 
backward direction of the polarization of the HCP and 
a rainbow phenomenon. Parts of these oscillatory struc- 
tures depend strongly on the ratio between the pulse du- 
ration of the exciting HCP and the classical Kepler period 
of the initially prepared Rydberg state. The amplitude 
and phase information which is contained in the semi- 
classical probability amplitudes might be useful for the 
reconstruction of quantum states of Rydberg electrons. 
A first indication of the usefulness of such an approach 
was given by the investigation of the spectra pertaining 
to radial wave packets. 
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FIG. 1. Classification of ionization spectra in the sud- 
den-impact approximation. The following semiclassical catas- 
trophes are observed in the respective areas of the plane 
spanned by Ap 2 /2|eo| and e^/|eo|: (I) forward and backward 
glory, (II) rainbow and backward glory, (III) two rainbows, 
forward and backward glory, (IV) no semiclassical catastro- 
phes. For a given value of Ap 2 /2|eo| the ionization proba- 
bility dPion/de^' is concentrated below the dashed line [in 
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FIG. 2. (a) Initial positions of classical trajectories with 
Ap 2 /2|e | = 6.25, e (/) /|eo| = 1-84, n = 50 according to 
Eq. (E7J). The numbers indicate final emission angles 0/ for 
specific initial positions. Angles are counted counterclock- 
wise from the z axis. The figure has to be continued into 
three dimensions by rotation around the z axis. The ro- 
man numbers relate the different trajectory classes to dia- 
gram (c). Axes are labelled in atomic units, (b) Final angle 
0/ as a function of the polar angle Qi of the corresponding 
initial position, (c) Classical weight = |P C !| 2 /|j4 c ;(xo)| 2 
(in a.u.) as a function of 0/ for the trajectory classes dis- 
tinguished in (a), (d) Angular distribution of the ionized 
electron ln{d 3 Pi 0rl /de'^' 1 df2[a.u.]} in the sudden-impact ap- 
proximation: quantum mechanical (full), uniform semiclassi- 
cal (dashed) and primitive semiclassical (dotted) result. 



FIG. 3. Same as Fig. | for e (/) /h| = 7.35. 
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FIG. 4. Angle- and energy-resolved ionization spec- 
tra obtained from radial Rydberg wave packets in the sud- 
den-impact approximation. (a) Outward travelling wave 
packets; the respective delay times between excitation and 
probing are indicated in units of the classical wave packet or- 
bit time T c i = 20 ps. (b) Incoming wave packets. Radial 
distributions of the wave packet at different times are shown 
in the insets. All relevant parameter values are given in the 
text. 
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HCP of duration T p = 



Same as Fig. |2j but now for the quasi-rectangular 
0.05T or 5 described in the text (de- 
flection function not shown), (a) The new trajectory classes 
which are denoted (IV) and (V) are separated by the "rain- 
bow point" G>. The initial points for (V) almost coincide with 
those for (III) (but the initial radial momenta are opposite). 
Indicated angles refer to (III), (c) Full curve: ionization prob- 
ability for the pulse described above. Dashed: sine-shaped 
pulse as described in the text. Dotted: result of the sud- 
den-impact approximation. 
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FIG. 6. Same as Fig. |(c) for e (/) /M = 3.68. The inset 
shows the spectra for the quasi-rectangular pulse described in 
the text (full curve) and for the sudden-impact approximation 
(dotted curve) at small angles. 



12 



